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Abstract 

We show that the phases of the 4-dimensional compact U(l) lattice gauge 
theory are unambiguously characterized by the topological properties of min- 
imal Dirac sheets as well as of monopole currents lines. We obtain the mini- 
mal sheets by a simulated-annealing procedure. Our results indicate that the 
equivalence classes of sheet structures are the physical relevant quantities and 
that intersections are not important. In conclusion we get a percolation-type 
view of the phases which holds beyond the particular boundary conditions 
used. 
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1. Introduction 



We investigate the 4-dimensional compact U(l) lattice gauge theory with Wilson 
action supplemented by a monopole term 

S = (5 Y, (! - cosG^) +\Y,\m p , x \ , 

p>u,x p,X 

where M p>x = e papv (Q) p ^ x+CT - Q^x)/^ and the physical flux Q^ x e [-it, it) is 
given by Q^x = Q^x + ^ n pv,x @- We consider periodic boundary conditions 
except for the discussion in Section 5 where we extend our study to a system with 
open boundary conditions. 

As is well known, the system has two phases, separated for A = by a first order 
phase transition. The strength of this transition decreases with increasing A until 
the transition ultimately becomes of second order ||. We have used this property 
to set up a very efficient algorithm 0, |3j] for Monte Carlo simulations in which A 
becomes a dynamical variable. The A dependence has also allowed us to study the 
dynamics of the transition in detail |§. 

We argue that the analysis of the topological structure of the configurations is a 
valuable tool for studying the phase structure, a tool which can provide much more 
detailed information than just the measurement of global observables. In the present 
work we use it by analyzing monopole currents and Dirac sheets. 

For individual loops the topological characterization is straightforward, but this 
becomes less trivial when loops are entangled in networks of monopole currents. 
Recently we have been able to produce a mathematically sound characterization of 
the latter M. This is not only necessary for the unambiguous identification of their 
physical features but also for the analysis of huge networks by computer. 

Here the analysis of configurations is extended to Dirac sheets and the fact that 
appropriately specified topological structures signal the phases is confirmed in more 
detail. In conclusion we point out a general principle for characterizing the phases, 
which does not rely on a particular choice of boundary conditions. 

2. Dual- lattice structures 

On the dual lattice the current J PjX = M p ^ x+P , defined over the links, obeys the 
conservation law Y^ P {Jp,x — Jp,x-p) = 0. Current lines are defined in terms of the 
current as follows: for J p x = there is no line on the link, for J p x = ±1 there is 
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one line, and for J PjX = ±2 there are two lines, in positive or negative direction, 
respectively. Networks of currents are connected sets of current lines. For a network 

— * 

N disconnected from the rest one can define a net current flow /, with components 

The Dirac string content of the plaquettes on the dual lattice is described in terms 
of a variable p pu>x = — \e papu n pv ^ x+p+(T , which satisfies the field equation ^2 a {p pc ,x ~ 
Ppa,x-a) — Jp,x- To every plaquette in the dual lattice we will associate no Dirac 
plaquette if p pa ,x = 0, one Dirac plaquette if p pcr ,x = ±1, with similar or opposite 
orientation according to the sign of p, and two Dirac plaquettes if p pa , x — ±2. 

Dirac sheets are formed by connecting Dirac plaquettes with a common edge and 
appropriate orientation, so that J PiX = except at the boundaries of the sheets. 
Dirac sheet structures are not gauge invariant. However, they belong to equivalence 
classes which cannot be deformed into each other by gauge transformations. For 
given boundaries there are topologically distinct possibilities for the sheet struc- 
tures. This is illustrated in Fig. 1 by a simple example in 2 dimensions, where 
a topologically nontrivial current network may be accompanied by a trivial or a 
nontrivial Dirac sheet. The equivalence classes of Dirac sheets thus carry more 
information than the related current networks. 

We represent the equivalence classes of the sheet structures by their members with 
minimal area, which are obtained by minimizing the number of Dirac plaquettes by 
a gauge transformation. This permits a unique specification and, in view of the 
large initial number of Dirac plaquettes, is necessary for a computational analysis. 

To address the connectedness of sheet structures, we construct minimal Dirac 
sheets by making first the connections where only two Dirac plaquettes meet at an 
edge with appropriate orientation. Then we consider the places where more than 
two Dirac plaquettes meet at an edge. If the sheets can be deformed (by a gauge 
transformation) in such a way that no more than two plaqettes meet we connect the 
respective pairs. Otherwise all must to be connected. 

In Fig. 2 we use some simple 3-dimensional examples to illustrate typical pos- 
sibilities for the intersection of Dirac sheets or for the meeting of more than two 
plaquettes. In Fig. 2a, deforming one of the sheets by a gauge transformation, the 
sheets can be separated. Consequently there is no connection. The same holds for 
Fig. 2b, where one must realize that common sites are not considered as a con- 
nection of sheets. In Fig. 2c, although one can identify two sheets, the structure 
cannot be separated by a gauge transformation and, therefore, is to be considered 
as connected. In Fig. 2d separation is obviously not possible. 
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Figure 1: Example of current network with two possibilities for a Dirac sheet on 
a 2-dimensional lattice with periodic boundary conditions. The shaded (unshaded) 
area represents a topologically non-trivial (trivial) sheet with given boundary, (b) 
shifted against (a) by L/2 in both directions. 
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3. Topological analysis 



The elements of the fundamental homotopy group 7Ti(X, 6) of a space X with 
base point b are equivalence classes of paths starting and ending at b which can be 
deformed continuously into each other. Its generators may be obtained embedding a 
sufficiently dense network N into X and performing suitable transformations which 
preserve homotopy. We observe that, if a given network N does not wrap around in 
all directions, then only the generators of a subgroup are produced. This provides 
an unambiguous characterization of networks. 

Choosing one vertex point of N to be the base point b and considering all paths 
which start and end at b, we note that a mapping which shrinks one edge to zero 
length preserves the homotopy of all of these paths. By a sequence of such mappings 
we can then shift all other vertices to b without changing the group content until 
we finally obtain a bouquet of paths which all start and end at b. 

Describing a path by a vector which is the sum of oriented steps along the path, 
for N with K vertices and K 1 edges the bouquet has K — K x — K + 1 loops, 
represented by vectors Sj, with i — 1, . . . , K and j-th components s^- = WijLj, where 
Lj is the lattice size. The bouquet matrix Wij must then be analyzed for the content 
of generators of 7Ti(T 4 , 6) = Z 4 . 

Current networks have the additional properties of path orientation and current 
conservation. The maps reducing the bouquet matrix Wij, therefore, in addition 
to homotopy have to respect current conservation. This leads to a modified Gauss 
elimination procedure, whereby the addition of a row to another requires its simul- 
taneous subtraction from a further row. In this way one arrives at the minimal form, 
with rows a±, . . ., a r , t, 0,. . ., where r < 4. The significance of this minimal matrix 

becomes obvious if one switches to the pair form with rows a±, —ai, . . ., a r , —a r , 
—»—»—» — » 
/, 0, . . ., 0, which exhibits the relation to the net current flow / explicitely. The 

number of independent pairs determines the number of nontrivial directions. 

For the topological analysis of minimal Dirac sheets we use the network of plaque- 
ttes as a sufficiently dense auxiliary network. The bouquet matrix then is obtained 
as described before. The reduction of the bouquet matrix is simpler here because 
the usual Gauss elimination applies, which gives the minimal form with rows a±, . . ., 

— # — # 

a r , 0,. . ., 0. Now the number of independent vectors corresponds to the number of 
nontrivial directions. 
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Figure 3: Probability -P ne t for the occurence of a nontrivial network in the hot and 
cold phase at the transition point as function of A for lattice sizes 8 4 (circles) andl6 4 
(crosses) . 

4. Numerical results 



Fig. 3 shows our results for the probability P net to find a current network which 
is nontrivial in four directions as function of A, with f3 chosen at the corresponding 
phase transition points || ffl. For larger L or negative A, where the peaks of the 
energy distribution related to the phases are well separated @, |J], P net is seen to 
be very close to 1 for the hot (confining) phase, and very close to for the cold 
(Coulomb) phase. We thus see that the two phases have an unambiguous topological 
characterization, provided by the existence of a nontrivial network in the hot phase 
and its absence in the cold phase. 

For A = the events with / ^ are rare; for lattice sizes L = 8 and L = 16 their 
percentage is less than 2% and 0.5%, respectively. For positive A they get slightly 
more frequent [§]. As is obvious from the pair representation of the bouquet matrix 
in Section 3, the net current flow / is not central for the topological characterization. 
In this context we should note that the quantity introduced in || to classify the 
networks, because of current conservation, is equal to /. 

For minimizing the number of Dirac plaquettes via a sequence of gauge transforma- 
tions we have used an annealing technique |J. Our procedure uses Metropolis sweeps 
based on a probability distribution P{rj) ~ exp(— arj), where rj = J2x,u>» \ n fiu,x\, 
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implemented through random local gauge transformations with angles maintained 
within [— it, it). In subsequent simulations we increase the parameter a in appro- 
priately chosen steps. First we determine r\ as function of a (which changes little 
with L and A) using tiny step sizes and large sweep numbers. Then we continue 
the procedure with 45 values of a (chosen in the interval [0.03,100] in such a way 
as to get approximately equal changes in if), performing up to 100 sweeps for each 
a. To be sure about the results, for each configuration we have generally performed 
three such runs as well as some additional runs with many more sweeps. Repeated 
applications of the method, with different random numbers, lead to consistent values 
for minimal numbers of Dirac plaquettes, with deviations smaller than 2%. Thus 
our procedure appears to produce a rather reliable determination of the absolute 
minimum. 

We found that the careful annealing described above is absolutely necessary. A 
minimization of the number of Dirac plaquettes based only on iterative gauge trans- 
formations gives strongly fluctuating results, which on the average are roughly fac- 
tors 2 and 1.3 larger in the cold and hot phase, respectively, than the minimal 
numbers obtained by our procedure. This is not only inappropriate in principle 
but also not sufficient for a proper computational identification of the topological 
structures. 

Before applying the annealing procedure the numbers of Dirac plaquettes are sim- 
ilar in both phases, but the numbers become markedly different after annealing. 
For example, for L = 8 and A = 0, in typical configurations with 8703 and 8793 
plaquettes in the hot and cold phase, respectively, our procedure reduced the num- 
bers to 306 and 846. The stronger reduction in the cold phase is mainly due to 
the disappearance of closed trivial structures. We also note that the procedure not 
only affects sheets separately, but involves the full set of them. For example, even 
pairs of oppositely oriented nontrivial sheets, first observed in ||, disappear with 
our annealing procedure. 

Our analysis of minimal Dirac sheet structures leads to results which agree with 
those derived from the study of current networks. The hot phase is characterized 
by the existence of a Dirac sheet which is nontrivial in all directions and the cold 
phase by its absence. To understand the significance of this result it is important 
to remember that, given a nontrivial current network, the sheets could be trivial or 
nontrivial (as illustrated by the example of Fig. 1). Thus our results show that in 
the hot phase the placement of the largest sheet is one which makes it topologically 
nontrivial. This points to the remarkable fact that the (equivalence classes of) Dirac 
sheet structures are the physically relevant objects. 
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Figure 4: Size of largest current network (circles, links) and largest minimal Dirac 
sheet (crosses, plaquettes) in the hot and cold phase at the transition point, as 
function of A on 8 4 lattice. 



In this context it should be remembered that some time ago Grosch et al. |7j were 
able to identify nontrivial Dirac sheets without boundaries as the structures respon- 
sible for the metastable states occurring under certain conditions in the cold phase. 
Here we have shown that nontrivial sheet structures with boundaries characterize 
the hot phase. 

To study the role of connectedness, we construct minimal Dirac sheets by making 
first the connections where only two Dirac plaquettes meet at an edge with appro- 
priate orientation. We note then that making either all or none of the connections 
where more than two plaquettes meet (with correct orientation) changes the num- 
bers of sheets and networks, and the sizes of the largest sheets, only by 10% or less 
and does not affect the topology. Looking into more detail we find that, for the most 
part, these changes are due to intersections of the type of Figs. 2a and 2b, where the 
sheets must actually be kept separated. Taking this into account the change reduces 
to 3% or less. We thus obtain the further remarkable result that intersections of 
minimal sheets are not important. 

In Fig. 4 we present the sizes of the largest current network and of the largest 
minimal Dirac sheet. Obviously the sheets give a more sensitive signal for the 
different phases. This may be explained by the fact that area instead of length 
enters. It is to be noted that the definition we have used for the connectedness of 
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Figure 5: Order parameters -P ne t and nmax/?Hot for A = as functions of j3 on 8 4 
lattice. 



sheets plays an important role in this result. If one considers as connected sheets 
which have just sites in common, then the signal provided by the size of the largest 
minimal sheets becomes comparable to the one given by the networks. 

Fig. 5 indicates that the probability P net of the occurrence of a non-trivial network, 
taking values very close 1 or for hot and cold phase respectively, may be a more 
sensitive order parameter than n max /n tot , the relative size of the largest network, 
advocated in 0. 
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5. Discussion 



The success of our topological order parameter confirms that the phase transition 
of the 4-dimensional compact U(l) gauge field theory concides with a percolation- 
type transition of the topological structures. Examples of theories with different 
behaviors have recently be discussed in 0. 

In order to see the general meaning of this results we note that a topologically 
nontrivial structure on a finite lattice with periodic boundary conditions, i.e. on the 
torus T 4 , corresponds to an infinite structure on an infinite lattice. Thus the hot 
phase is more generally characterized by the existence of an infinite current network 
or Dirac sheet and the cold phase by its absence, where on finite lattices "infinite" 
is to be defined in accordance with the particular boundary conditions which are 
used. 

To test this characterization we also performed simulations with open boundary 
conditions [^Oj] . In the corresponding analysis the prescription "nontrivial in all 



directions" is replaced by "touching the boundary in all directions". In this way 
we again get an unambiguous signal for the phases. The comparison of the order 
parameters looks very similar to Fig. 5, the advantage of our order parameter being 
even more pronounced. Because of finite-size effects, however, the width of transition 
region becomes larger, for L = 16, by a factor of approximately 100 . 

In order to obtain information about the order of the phase transition the most 
immediate thing to do is to look whether there is an energy gap for sufficiently large 
L. With periodic boundaries L = 8 is sufficient for seeing a gap whereas L = 4 
is not. With open boundaries a gap is not yet seen for L = 16. The width of 
the peak is relatively narrow (comparable to the width of the peaks with L = 8 
in the periodic-boundary case) [fTCfl . It appears to us that lattices well larger than 
L = 16 will be needed in order to see a possible gap in systems with open boundary 
conditions. 



In [[yj the simulations have been performed on the surface of a 5 dimensional 
cube which is homeomorphic to the sphere S 4 . Considering lattices up to L = 10 
the authors observe only one peak. Its width for L = 10 is comparable to those of 
the peaks on T 4 for L = 8. The extension of the transition region, which could also 
provide information about finite-size effects, has not been reported. 

In [[12] simulations for lattices up to L = 16 have been done with fixed boundaries 
(i.e. setting the group elements to 1 at the boundaries), which makes the lattice 
again homeomorphic to S 4 . With these boundary conditions the transition region 
was observed to be wider by a factor of about 10 3 , which reflects the huge finite-size 
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effects caused by the strong inhomogeneity of the system. Because of low statistics 
the results in this work are of semiqualitative nature and the energy distribution is 
not given. 

It seems to us that in the two geometries homeomorphic to S 4 mentioned above 
larger lattices would also be necessary for deciding about the existence of a gap. 
Further, in these geometries the analysis of configurations of monopole currents and 
Dirac sheets would be of great interest, too. Appropriate definitions of "infinite" 
on finite lattices should of course be given. It should be realized that these geome- 
tries are different from a physical point of view. Fixed boundary conditions can be 
expected to lead to the same (flat-space) limit as periodic or open boundary condi- 
tions, while for the surface of a 5-dimensional hypercube some curvature (depending 
on the details of the limit) could persist. 

A further type of boundary condition, proposed in fl3|| , consists in the suppression 
of (all) monopoles at the boundaries. This is done by using an action leading to 
the suppression of monopoles |T3J] at the boundaries and the usual Wilson action 
elsewhere. Though the definition of "infinite" gets somewhat dilute in that case we 
expect it to fit into our picture as well. The existence of a gap with this type of 
boundary condition remains to be confirmed. 
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